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13,  ABSTRACT  (Mammurmoo  words) 

Much  of  the  extensive  literature  on  shell  theory  centers  around  the  computation  of  dispersion  curves,  usually  thought 
of  as  the  roots  of  a  determinant  D  which  relates  the  angular  frequency  w  to  the  wave  numbers  for  free  waves  traveling  on 
the  shell.  For  a  cylindrical  shell,  the  free  wave  paths  are  actually  helical  curves  on  the  surface  of  the  shell,  having  both  a 
circumferential  and  an  axial  wave  number.  In  this  case,  it  is  useful  to  indicate  explicitly  that  the  dispersion  relation  is  a 
function  of  three  variables,  i.e.,  m,  n)  =  0,  where  m  is  the  number  of  half  wavelengths  of  a  wave  traveling  in  the  axial 
direction,  and  n  is  the  number  of  fuU  wavelengths  of  a  wave  traveling  around  the  shell.  We  shall  refer  to  this  function  D(o), 
m,  n)  as  the  dispersion  volume  density,  since  it  is  a  function  of  three  variables.  Normally,  one  sees  various  two-dimensional 
displays  of  the  dispersion  relation,  which  are  created  by  plotting  any  two  of  these  variables  against  each  other,  while  hold¬ 
ing  the  third  variable  fixed. 


93“ 1 1429 

iilllllillliilllllillllilit 


Published  in  Proceedings,  14th  Internationa]  Congress  on  Acoustics,  Vol.  1,  3-10  Sep  1992,  paper  B8-2. 


14  SUBJECT  TERMS 


target  .strength 


martsam 


15  NUMBlBOPPAGCS 


16  PRICE  CODE 


17.  SECURITY  CLASSJPICATtON 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20  UMlTATiON  or-  ABSTRACT 


SAME  AS  REPORT 


NSN  7540-01 -ZeO-SSOO 


Standard  form  298  (BACK) 

UNCLASSIFIED 


B8-2 


i  nr.  KKKiriKS  I  TALCJ  l.A!  {UN  ANIi  DISJM.AV  OK  DISKKJISION 
rrH\  {-S  A  U{[N  I  IMJfifCAI.  SHKI.{.  {MMKItNKD  IN  A 

Finn 

Harry  A  SCHFNTK 

lU)r  ari<i  K  D!\i;^ion.  Naval  ('omnian<l  and  Ocfan  Surveillance 

Onter,  San  <^A  92152-5000 

Much  of  ilip  extensive  Hlerature  on  sh*'!!  theory  centers  around  the 
'f  dispersion  curves,  usually  *hou^hl  of  as  the  roots  of 
a  deterr'iir'''’an(  H  «'hKh  re)a‘‘'s  ^■»'qiiency  u»  to  the  wave 

nuriihers  for  fr-e  waves  traveling  on  th*-  slie}!  For  a  cylindrical  shell, 
the  free  wa\e  paths  are  actually  helical  curves  on  the  surface  of  the 
shell,  having  both  a  circumferential  and  an  axial  wave  number.  In  this 
case,  it  is  u.‘ieful  to  indicate  explscUly  that  the  dispersion  relation  is  a 
function  of  three  variables,  i.e,,  0(w,m,ti)  =  0.  where  rn  is  the  number 
of  half  wavelengths  of  a  wave  traveling  in  the  axial  direction,  and  n  is 
the  number  of  full  wavelengths  of  a  wave  travelir  *  around  the  shell-  We 
shall  refer  to  this  function  £)(«*/, rn.n)  as  the  dispersion  volume  densifjr. 
since  it  is  a  function  of  three  variables  Normally,  one  sees  various  two- 
dimensional  displays  of  the  dispersion  relation,  which  are  created  by 
plotting  any  two  of  these  variables  against  each  other,  while  holding  the 
third  variable  fixed,  as  illustrated  m  Fig  1 

Dispersion  relations  for  a  fluid-loaded  cylindrical  shell  have  been 
developed  previously  by  many  authors  .  However,  as  pointed  out  in  a 
lengthy  article  by  Scott  (1],  some  of  these  pre  viously  published  works 
(e.g.,  ^2]  and  made  the  mistake  of  searching  for  roots  of  the  real 

part  of  the  dispersion  equations,  which  led  to  the  computations  of  spu¬ 
rious  non-physical  branches  of  the  dispersion  curves.  Although  Scott 
has  convincingly  made  this  point,  his  actual  computation  procedure 
is  cumbersome.  In  fact,  it  is  difficult  to  determine  the  complex  roots 
of  such  dispersion  relationships  by  any  root-following  method  la  the 
present  work,  the  complications  of  calculating  and  following  these  com¬ 
plex  roots  have  been  simplified  by  using  a  regular  r.jmputational  grid 
for  each  parameter  of  the  dispersion  function,  and  th'»n  making  the  root 
loci  evident  by  tht  appropriate  use  of  graphics  to  display  what  we  term 
the  dispersion  tv/umr  denstty. 

The  purpose  of  this  paper  is  to  describe  and  illustrate  an  efficient 
procedure  for  calculating  and  displaying  the  dispersion  curves.  We  will 
base  the  results  on  Flugge  shell  theory,  which  ha<:  some  advantages  for 
air-backed  thin  steel  shells  immersed  in  water.  Only  the  critical  steps 
and  results  in  the  derivation  of  the  dispersion  relationships  will  be  given 
here 


The  displacements  and  forces  are  expanded  in  a  combined  Fourier 
senes  and  transform  to  allow  for  completely  general  motions  due  to 
asymmetric  for^-ing  functions.  With  a  time-harmonic  factor  of  <’"*  sup¬ 
pressed,  the  displacements  of  the  rnid-surface  of  a  straight¬ 

sided  cylindrical  shell  of  mean  radius  a  and  thickness  h  are  a  function 
of  i  and  only.  In  general,  the  forcing  function  F  will  also  have  three 
degrees  of  freedom,  but  F  will  be  identically  zero  for  a  shell  in  eacao, 
and  F*  =  F*  ~  0  when  the  only  force  causing  or  resisting  motion  of 
the  shell  is  p(a.  d,  z),  the  total  (incident  plus  scattered)  pressure  on  the 
shell  due  to  an  incident  acoustic  wave  The  shell  material  will  be  rhar- 
aclenzed  by  three  parameters.  Poisson's  ratio  u,  Young's  nwduhis  E. 
and  the  density  p,;  the  fluid  by  a  density  pj  and  a  speed  of  sound  C/ 
Other  (non-dimensional)  parameters  describing  the  solution  arc 


CircumfercnliaJ  wave  number,  n 
(c)  FYequency  Q  vs  ha/rmunr  n  for  fixed 
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We  assume  that  u'(^,b)  can  be  i-xpanded  m  a  Fourier  senes  in  d>. 


(aj  Conventional  dispersion  plot,  <;  vs  Q  far  fixed  i 


tion  / 


Fig  1  (jits  along  the  principal  axes  of  Ihr  <ji«pcr.«ion  »<>/vrne  dcn.irty  “ 

for  a  thin  rylindricaJ  shell  in  water  vX^df 
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Avfill  j'id/Of 


,i\-\  ^  ^0 


(1) 


and  that  the  axial  variation  of  »»(*)  expreaaed  as  a  Fourier 

transformation 


<u)  =  r 

J-on 


where  the  transform  variable  (  is  a  non-dimensional  axial  wavenumber. 
Similar  expansions  are  assumed  for  u^(^,s)  and  ti'(^,s).  Recause  of 
the  orthog<ma3ity  properties  of  the  Fourier  expansions,  the  application 
of  the  operators  to  the  expansions  for  the  displacement  results  in  an  in¬ 
finite  set  of  matrix  equations  for  each  one  for  each  harmonic  order  n. 
The  differential  operators  all  become  matrix  operators,  resulting  in  a  set 
of  coupled  linear  equations.  We  also  represent  the  acoustic  pressure  in 
a  double  Fourier  expansion  which  allows  us  to  express  the  forcing  term 
as  a  function  of  U'(0>  ^He  double  Fourier  expansion  of  the  radial  dis¬ 
placement.  This  permits  us  to  move  the  forcing  term  in  the  transformed 
equations  over  to  the  left-hand  side  and  write  a  completely  honnogeneous 
matrix  equation  for  the  Fourier  components  of  shell  displacement; 


Koo  Koi  Kqj 

0i«)  f  0  ] 

KU=  /f.o  Ku  Ki,  VHO  =  0  . 

10. 

.  /fjo  /fjl  /fjj  ^ 


The  components  of  K  in  equation  (3)  are  given  by 


Koo  =  -<»  -  [^n’ 


Koi  s=  Kio  =  -n<- 


Koi  =  Kao  =  j ^ 


11  = - 5 — s  - s — s 


Kia  =  Kai 


,  .  Pl^^f 


in  which  the  additional  terms  to  account  for  the  resistance  and 
,  reactance  components  of  the  fluid  loading  are 


»•««)=«/>/«  Im 


»HP(V)J' 


and  y  «  \/(to)*  -1;*  for  frequencies  such  that  ka  ^  ^  {k  Z  *1)- 
*Whsn  ka<({k  <  k,t  lh<  lo-eslled  tvsneietnt  region), 


where  r  ss  -  (ifcfl)*- 

Kquation  (3)  can  have  a  noti-lrivial  solution  whenever  its  determi¬ 
nant  vanishes,  which  can  occur  for  particular  combinations  of  frequency 
n,  axial  wavenumber  C>  snd  circumferential  wavenumber  n  when  there  is 
no  fluid  loading.  These  combinations  represent  the  dispersion  relations 
for  this  shell  in  vacuo.  Fluid  loading  introduces  (acoustic)  damping  or 
loss.  Even  in  this  case,  the  absolute  value  of  the  determinant  can  be¬ 
come  very  small,  and  the  resultant  forced  solution  can  be  very  large  for 
particular  conibinalions  of  0,  and  n.  Thus  minima  of  the  absolute 
value  of  the  determinant  of  K  will  define  the  dispersion  relations  for  the 
fluid-loaded  ease. 

The  simplified  computational  procedure  is  to  compute  the  coeffi¬ 
cients  of  the  matrix  K  and  then  evaluate  the  determinant  of  K  at  regular 
increments  in  the  principal  variables  which  are  the  non-dimensional  fre¬ 
quency  n,  axial  wavenumber  (,  and  circuniferenliai  wavenumber  n.  The 
power  of  the  method  lies  in  the  fact  that  modern  graphical  workstations 
with  appropriate  visualization  software  can  not  only  perform  these  cal¬ 
culations  but  can  quickly  display  the  results  in  a  manner  which  allows 
one  to  understand  the  dispersion  relations  in  a  global  sense  (over  the 
entire  meaningful  range  of  0,  and  n.  The  actual  program  code  was 
written  in  an  efficient  high  level  mathematical  language. 

Once  the  dispersion  volume  density  is  calculated  for  an  appropriate 
range  and  density  of  0,  C*  end  n,  the  results  can  be  displayed  in  a  variety 
of  ways.  One  of  the  simplest  and  moat  effective  means  it  compute  the 
logarithm  of  the  reciprocal  oi  the  absolute  value  of  D,  and  then  create 
a  smooth  interpolated  color  image  of  the  results.  If  the  color  table  is 
properly  chosen,  Ibis  has  the  effect  of  showing  the  root  loci  as  lines  whose 
brightness  indicates  the  nearness  to  zero  and  whose  width  is  indicative 
of  width  of  the  null  and  of  the  resolution  in  the  appropriate  parameters 
(0,  C>  ^)‘  If  calculated  in  this  way,  it  is  rJso  possible  to  animate  the 
display  of  these  color  images,  which  enhances  the  understanding  of  bow 
the  location  of  the  modes  change  with  changes  in  various  parameters. 
All  of  these  ideas  will  be  illustrated  in  the  presenfefisn  of  ibis  paper, 
but  cannot  be  adequately  shown  in  this  brief  black  and  white  hardcopy. 

In  summary,  the  present  work  confirms  the  earlier  work  of  Scott 
as  to  the  correct  dispersion  loci,  but  provides  a  simpler  and  more  easily 
interpreted  method  of  calculating  and  displaying  the  dupersien  volume 
ienoitf  over  the  entire  relevant  range  of  frequency  and  wavenumbers. 

This  work  wss  supported  by  the  Naval  Research  Laboratory  and 
the  Office  of  Naval  Technology.  The  suggestions  and  critique  of  the 
mathemi.lics  and  computalicmal  procedures  by  George  W.  Denlhien  is 
grstefully  acknowledged. 
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In  the  original  submission  of  this  paper,  a  factor  of  was  omitted  from  several 
equations.  Equations  (6)-(8)  of  Paper  B8-2  should  read  as  follows; 


mim 


Ko2  =  K20  =  -iCi'  -  (6) 

Ku  =  +  {kpaf  -  (7) 

K,2  =  K21  =  -in  -  (8) 


